AD-A222  839 


fifLQ 


DTIO  FiLE  CORY 


The  Survival  Probability  Function  Of  A  Target  Moving  Along  A 
Straight  Line  In  A  Random  Field  Of  Obscuring  Elements 

by 
•  v 

S.  Zacks  and  M.  Yadin 

May  15,  1990 

Technical  Report  No.  1 

Prepared  Under  Contract  DAAL03-89-K-Q129  with  the  U.S.  Army  Research  Office 

S.  Zacks,  Principal  Investigator 

Center  For  Statistics,  Quality  Control  and  Design 
State  University  of  New  York 
University  Center  at  Binghamton 
Binghamton,  NY  13901 


-'V 


&  JUN  1 5  1390 $  * 


& 


f'.k’T"  STATEMENT  A 

Apprised  for  public  relaas* 
Wamour.cn  Unlnxuted 


I 


The  Survival  Probability  Function  Of  A  Target  Moving  Along  A 
Straight  Line  In  A  Random  Field  Of  Obscuring  Elements 


by 


S.  Zacks  and  M.  Yadin 
State  University  of  New  York  at  Binghamton 

and 

Technion,  Israel  Institute  of  Technology 


Abstract 

A  target  is  moving  along  a  straight  line  path.  Random  portions  of  the  path  might 
be  invisible  to  the  hunter  (in  shadow).  Shooting  trials  are  conducted  only  along  the 
visible  segments  of  the  path.  An  algorithm  for  the  numerical  determination  of  the  survival 
probability  of  the  target  is  developed.  This  algorithm  is  based  on  the  distribution  of 
shadow  length  which  is  also  developed. 


Key  Words :  Lines  of  sight;  visibility  probabilities;  distributions  of  shadows;  survival  prob¬ 
ability. 

Acknowledgement:  Research  supported  by  the  U.S.  Army  Research  Office,  Contract  DAAL03- 
89-K-0129  with  the  State  University  of  New  York  at  Binghamton. 


1.  Introduction. 

„ _ T>  r 

Tho-pfe-iLiTt  study  is  focused  on  the  problem  of  determining  the  survival  probability  of 
a  moving  target,  which  is  under  attack  by  a  hunter.  The  target  (vehicle,  tank,  etc.)  is 
moving  along  a  straight  line  path,  which  is  partially  obscured- from  the  hunter  by  randomly 
distributed  objects  (trees,  clouds,  terrain  objects,  etc.).  The  target  can  be  destroyed  by  the 

hunter  only  along  the  visible  segments  of  the  path.  Visibility  contact  between  the  hunter 

f GlM-  b  O 

and  the  target  is  needed  for  ^  time  units  for  a  shooting  trial  to  occ  ar.  In  any  given  shooting 

trial  the  probability  that  the  target  is  destroyed  is  fixed.  If  the  target  survives  a  shooting 

trial,  another  identical  trial  may  be  attempted  if  continuous  visibility  for  tq  time  units  is 

possible.  If  the  target  enters  an  obscured  segment  of  the  path,  the  shooting  trials  terminate, 

until  visibility  contact  is  reestablished.  Under  the  above  assumptions,  if  the  target  has  to 

cross  a  visible  segment  of  length  L ,  its  survival  probability  can  be  approximated  by  the 

t'n  r  _ * 

negative  exponential  function  exp{—  qL)  ,  for  suitably  chosen  constant  g,  0  <  q  <&(rrThe 

problem  is  that  the  number  of  visible  segments  on  the  moving  path,  between  two  specified 

1-  SVfc  O- 

points  pt  and  pf.  and  their  lengths  are  random  variables,  whose  distributions  depend 
on  the  characteristics  of  the  random  field, 

C  '  Tfrr  ~  u'T  ~  , 

TJt*  UF"Cjprii:  study  is  based  on  the  model  of  a  random  Poisson  field  of  obscuring  elements. 
This  model  is  presented  in  Section  2.  Under  the  assumptions  of  this  model,  it  is  relatively 
simple  to  derive  the  conditional  distribution  of  the  length  of  a  visible  segment  on  the  right 
hand  side  (r.h.s.)  of  a  point,  Px ,  on  the  path,  given  that  the  point  Px  is  visible.  This 
distribution  is  given  in  Section  3.1.  On  the  other  hand,  it  is  more  complicated  to  determine 
the  distribution  of  the  length  of  a  segment  which  is  obscured  (in  shadow).  In  Section  3.2 
we  present  the  methodology  for  determining  the  distribution  of  the  length  of  shadows. 
This  methodology  is  based  on  a  theory  given  by  Chernoff  and  Daly  (1957).  For  a  given 
point  Px  on  the  moving  path,  Chernoff  and  Daly  (C-D)  define  the  functional  T(x ),  which 
is  the  right  hand  limit  of  the  shadow  to  the  r.h.s.  of  Px ,  cast  by  obscuring  elements  in 
the  field  which  intersect  the  ray  Rx  from  the  origin  O  through  Pz  .  Employing  functions 


K±(x ,  y  '■  which  are  defined  in  Section  2  and  derived  explicitly  for  the  standard-uniform 
case  in  the  Appendix,  we  express  the  cumulative  distribution  function  (c.d.f.)  of  T(x) 
explicitly.  The  right  hand  limit  of  a  shadow  to  the  r.h.s.  of  Pr,  is  U(x )  =  lim  Tn(x), 

n — "oo 

where  Tn+1(x )  =  T(T”(x)),  for  n  =  0, 1, ,  T°(x)  =  x.  The  relationship  between 
the  c.d.f.  of  Tn+1(x)  to  that  of  Tn(x),  n  =  0,1,*"  is  discussed  in  Section  3.2.  The 
distribution  of  U(x)  is  obtained  as  a  limit  of  that  of  Tn(x).  From  the  distribution  of  U(x) 
we  obtain  the  conditional  distribution  of  the  right  end  limit  of  a  shadow  to  the  r.h.s.  of 
Px ,  given  that  Px  is  the  first  point  in  the  shadow. 

In  Section  4  we  employ  the  results  of  Section  3  to  approximate  the  survival  probability 
function  5(x,y)  along  the  moving  path  between  the  points  P*  and  Pa ,  x  <  y.  The 
function  S(x,y)  is  given  by  the  integral  equation 

(1.1)  S(x,y)  =  A(x,y)  +  J  B(x,w)S(w,y)dw 

where  A(x,y)  and  B(x,y)  are  defined  in  terms  of  the  distributions  of  the  lengths  of 
visible  and  non- visible  random  segments,  as  shown  in  Section  4.  An  algorithm  for  the 
discrete  approximation  of  the  solution  of  (1.1)  is  given  in  Section  5.  Numerical  solutions 
based  on  this  algorithm  are  provided  there  too.  A  Quick  Basic  program  (version  4.5)  for 
computations  can  be  obtained  upon  request. 

In  a  previous  Technical  Report  [7]  we  approximated  the  survival  probabilities  by  deriving 
lower  and  upper  bounds  to  the  distribution  of  the  number  of  shooting  trials,  N,  along 
the  path.  The  present  study  provides  the  method  of  computing  the  survival  probability 
function  5(x,y),  which  is  required  for  various  applications.  W‘*h  the  new  algorithms  for 
determining  distributions  of  shadows  and  survival  functions  we  i  tackle  problems  like 
the  Hunter-Escort  problem,  which  will  be  discussed  in  another  paper. 
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2.  The  Random  Field  Model  And  The  Determination  Of  Visibility  Probabilities. 

In  the  present  paper  we  consider  a  two  dimensional  physical  model.  A  generalization  to 
a  three  dimensional  model  can  be  done  in  a  similar  fashion  to  that  of  Yadin  and  Zacks  [5]. 
The  moving  path  of  the  target  is  a  straight  line  C.  The  Hunter  is  located  at  a  point  O  (the 
origin),  at  distance  r  from  C.  Let  U  and  W  be  two  straight  lines  parallel  to  C ,  located 
between  O  and  C,  at  distances  u  and  w  from  O,  respectively;  0  <  u  <  w  <  r.  The 
obscuring  objects  are  modeled  by  a  countable  number  of  disks  of  random  size,  which  are 
centered  at  random  points  in  the  strip  5,  bounded  by  U  and  W.  We  consider  a  cartesian 
coordinate  system  in  which  the  y-axis  is  a  straight  line  through  O,  perpendicular  to  C, 
which  intersects  ZV,  W  and  C  at  the  point  (0,  u),  (0,  w)  and  (0,r),  respectively.  A  point 
Pr  on  C  has  coordinates  (x,  r) . 

A  random  disk  is  represented  by  the  random  vector  (X,Y,Z),  where  (X.Y)  are  the 
random  coordinates  of  the  center  of  the  disk,  and  Z  is  its  random  radius.  Without  loss  of 
generality,  assume  that  the  sample  space  of  (X,Y,Z)  is  S*  =  S  x  [a,  6],  where  0  <  a  < 
b  <  o o.  Let  B*  be  the  Borel  a  -field  on  S* .  Let  {(X;,  Fj,  Z,) ,  i  =  1,2,  •••}  represent  a 
sequence  of  countable  random  disks  measurable  w.r.t.  the  same  space  (S* ,B* ,  P) .  It  is 
assumed  that  the  random  vectors  are  independent  and  identically  distributed  (i.i.d.),  and 
have  a  common  distribution  H(x ,  y,  z) .  Let  F(z  j  x,  y)  denote  the  conditional  c.d.f.  of  the 
radius  Z,  given  the  center  (X,  Y)  is  at  (x,y).  Let  h(x,y)  be  the  joint  p.d.f.  of  (X,  Y), 
such  that  h(x,y)  =  0  for  all  (x,  y)  £  S.  We  further  assume  that  the  probability  that  a 
random  disk  intersects  either  O  or  C  is  zero.  Let  B  be  any  Borel  set  in  B* .  Let  N{B} 
designate  the  number  of  random  disks  with  coordinates  in  B . 

If  {B i,--  -  ,Bm}  is  any  finite  partition  of  5*,  m  =  1,2, •••  ,  it  is  assumed  that  the 
random  variables  7V{5,},  i  =  l,--  -  ,m  are  independent,  having  Poisson  distributions 
with  expected  values 

(2.1)  = 

Bi 
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0  <  A  <  oo.  Such  a  random  field  is  called  a  Poisson  field.  The  Poisson  field  is  called 
standard-uniform  if  dH(x,y,  z)  =  hlc(x,y)f(z)dxdydz,  where  0  <  h  <  oo ,  C  is  a  subset 
of  S  which  represents  the  field  of  view  of  the  Hunter,  and  Ic(x,  y)  is  the  indicator  function 
of  C.  A  point  Pz  on  C  is  said  to  be  visible  from  O ,  if  the  ray  Rz  from  O  through  Pz  is  not 
intersected  by  random  disks.  In  a  similar  manner  we  can  define  the  notion  of  simultaneous 
visibility  of  several  points  on  C.  In  our  previous  papers  [2,3,4]  we  have  introduced  the 
functions  K+(x,t)  and  K-(x,  t)  for  0  <  t  <  oo;  where  \K±(x,  t)  is  the  expected  number 
of  disks  centered  in  S  between  the  rays  Rz  and  Rz±t,  which  do  not  intersect  Rz  .  Explicit 
formulae  for  K±(x,t),  for  the  standard-uniform  case,  with  a  uniform  distribution  of  radii 
on  [a,  6],  0  <  a  <  b.  is  given  in  the  appendix. 

Let  [L,  U]  be  an  interval  of  the  x  -coordinates  of  the  point  on  C  belonging  to  a  segment 
of  interest.  Let  L*  <  L  and  U*  >  U  be  properly  chosen,  and  C*  the  set  in  S  (trapez) 
between  the  rays  R  and  Ry.  .  One  can  verify  that  the  probability  that  Pz  is  visible, 
for  some  L  <  x  <  U ,  is 

(2.2)  0(x)  =  exp{-[fi{C*}  -  -  L*)  -  \K+(x,U*  -*)]}. 

For  a  formula  of  the  simultaneous  visibility  of  n  points  in  [ L ,  17] ,  see  Yadin  and  Zacks  [4]. 


3.  Distributions  Of  Length  Of  Visible  And  Of  Shadowed  Segments. 

3.1.  Distributions  Of  The  Length  Of  Visible  Segments. 

In  the  present  section  we  derive  a  formula  for  the  conditional  c.d.f.  of  the  length  of  a 
visible  segment  to  the  r.h.s.  of  Pz ,  given  that  Pz  is  visible. 

Let  I(x)  be  an  indicator  function  which  assumes  the  value  1  if  Pz  is  visible,  and  the 
value  zero  otherwise. 

Let  L(x)  be  the  length  of  the  visible  segment  of  C  to  the  r.h.s.  of  Pz ,  i.e., 

(3.1)  L(x)  =  inf{y  :  y  >  x,  JJ  I{u)  =  1}  -  x. 

x<u<y 
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We  derive  here  the  formula  for 

V{1  |  x)  =  P{L(x )  <  l  |  I(x)  =  1}. 

(3.2) 

=  1  -P{L(x)  >  l  |  I(x)  =  1}. 

Let  C*  be  the  set  of  (x,y)  points  in  <S,  which  was  defined  in  the  previous  section.  We 
derive  the  formula  of  V(l  |  x),  for  L  <  x  <  U ,  and  0  <  l  <  u  —  x . 

Let  C_(x)  be  the  set  bounded  by  U ,  W  and  the  rays  Rf  and  Rx  .  Let  C(x,l)  be 
the  set  bounded  by  U ,  W  and  the  rays  Rz  ,  Rx+j;  and  C+(l  +  x)  the  set  bounded  by  U . 
W,  R;+r  and  Ry  .  Notice  that  C*  =  C-(x )  U  C(x,l )  U  C+(l  +  x).  As  before,  we  denote 
by  /u{C }  the  expected  number  of  disks  having  centers  at  the  set  C,  as  given  by  (2.1). 
Accordingly, 

(3.3) 

P{L{x)  >  /,  I(x)  =  1}  =  exp{-[/i{C_(x)}  -  AA'_(x,x  -  L *)]  -  n{C(x,  /)} 

-  [n{C+(l  +  x)}  -  A K+(l  +  x,  U*  -  l  -  x)]} 

=  exp{— /x{C*}  +  A[A'_(x,  x  —  L*)  +  I\+(l  +  x,U*  —  l  —  x)]}. 
Dividing  (3.3)  by  (2.2)  we  obtain 

(3.4)  P{L{x)  >  l  |  I(x)  =  1}  =  exp{-A[A+(x,  U*  -  x)  -  A'+(/  +  x,U*  -  l  -  x)]}. 

3.2.  The  Distribution  of  Shadow  Length. 

We  have  denoted  by  U(x)  the  right  hand  limit  of  the  shadow  on  C  to  the  r.h.s.  of  Pr . 
Let  D(u  |  x)  denote  the  conditional  c.d.f.  of  U(x),  given  that  the  shadow  starts  at  P* . 

Consider  the  rays  Rr  and  Ry  for  y  >  x.  Let  N(x,y )  denote  the  number  of  disks 
centered  in  S ,  which  intersect  both  Rr  and  Ry.  Define  the  functional 

(3.5)  T(x)  =  sup{y  :  N(x,y)  >  1}. 

Furthermore,  let  T,+  1(x)  =  T(T‘(x)),  i  =  0, 1,  •  •  •  where  T°(x)  =  x.  Obviously,  T,+1(x)  > 
T'(x) ,  for  all  i  >  0,  and  therefore  U{x)  =  limT‘(x).  U{x)-x  is  the  length  of  the  shadow 
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Co  the  r.h.s.  of  Px  .  We  derive  first  the  c.d.f.  of  T(x).  Clearly,  { T(x )  >  t}  =  {.V(x,  t)  >  1}  . 
Thus, 

(3.6)  P{T(x)  <  t)  =  P{N(x,t )  =  0}  =  exp{—/u(x,t)}, 

where  =  E{N(x,  t)} .  Furthermore, 

(3.7)  t)  =  fi{C*}  —  XK+(x,  U*  —  x)  —  AA'_(t,  t  —  L*)  +  \K+(x,t  —  x)  4-  A t  —  tn 


where  t  is  the  coordinate  of  the  bisector  between  Rj,  and  R* .  Notice  that,  since  A’+(x,  0)  = 
A'_(x,0)  =  0  for  all  x, 


(3.S) 


p(x,  x)  =  limp(x,  t) 

tlx 

=  /i{C*}  —  AA'+(x,  u*  —  x)  —  AA‘_(x,  x  —  L*). 


Hence, 


(3.9)  li  mP{T(x)<t}  =  if>(x), 

t  Lx 

which  is  the  probability  that  Pz  is  visible.  Thus,  the  c.d.f.  of  T(x),  H(t-,x)  is  zero  for 
t  <  x,  it  has  a  jump  point  at  x,  H(x]x)  =  ip(x),  and  is  absolutely  continuous  at  t  >  x. 
This  property  is  inherited  by  the  c.d.f.  of  T”(x),  i7n(t;x).  We  provide  now  the  recursive 
relationship  between  Hn(t]x )  and  Hn-\(t;x).  Introduce  the  bivariate  distribution 

f2;*)  =  P{Tn~l{x)  <  h,Tn(x)  <  t2}. 
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Indeed,  given  that  {Tn  2(x)  —  z,Tn  *(x)  =  y},  {Tn{x)  >  f}  if  and  only  if,  there  exists 
at  least  one  disk  which  intersects  and  Rt,  but  does  not  intersect  R-.  Hence, 


/•n  yt  i 

(3.12)  Gn(t1,t2;x)  =  J  J  exp{-[n{u,t2)  -  /j(z,t2)}}dG„-i(z,u;x). 

These  bivariate  c.d.f.  can  be  determined  recursively,  starting  with  Gi(t\ ,  t2;  x)  =  H(t2\x) 
for  all  1 1  <  t2 .  Moreover, 

(3.13)  G2(tut2;x)  =  J  1  Qf  e^'^dz'j  dH{u-x). 

Finally,  since  Hn+i(t;x )  <  Hn{t\x)  for  each  t  >  x  and  all  n  =  1.2,  the  c.d.f.  of  U(x) 


is 


:  3.14) 


P{U(x)  <  t}  =  lim  Hn(t;x). 


Thus.  P{U(x)  <  t}  =  0  for  all  t  <  x,  and  lim  P{U(x)  <  t}  =  ?/>(x).  The  conditional 

Ux 

c.d.f.  of  U(x).  given  {/(x)  —  0}  is 


(  P{U(x )  <  <}  - 


(3.15) 


P{U(x)  <  t  |  I{x)  =  0}  =  { 


,  for  t  >  x 


1  - 

V  0,  for  t  <  x. 

We  are  interested,  however,  in  the  conditional  c.d.f.  D(u  |  x),  where  Px  is  the  first  point 
(the  left  hand  limit)  of  the  random  segment  in  shadow. 

Simple  geometric  considerations  yield  that  the  length  of  a  random  shadow  cast  by  a 
single  disk,  having  left  hand  limit  at  P* ,  with  center  on  a  line  parallel  to  U  at  distance  h 
from  O ,  and  disk  radius  Z ,  is 


(3.16) 


U(x,  h,Z)  =  r  tan  2  sin 


-l 


+  tan' 


■'(?) 


-  x, 


where  (xc,  h)  are  the  coordinates  of  the  center  of  the  disk,  with 

1/2 


(3.17) 


xc  =  —  h  Z  [  1  + 
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Thus,  if  a  <  Z  <  b  w.p.  1,  the  minimal  length  of  shadow  starting  at  x  is 
(3. IS) 


Finally,  since  a  shadow  starting  at  Pr  ends  at  point  U{x)  >  um(x)  +  x , 


(3.19) 


^  N  P{U(x)  <u}~  P{U(x)  <  um(x)  +  x} 

D(u]x)  = - ~f{gW<i-wV.) - 


for  u  >  x. 


In  Section  5  we  provide  computing  algorithms  for  the  numerical  determination  of  the 


xd.f.’s  V(l  |  x)  and  D(u  |  x),  and  illustrate  them  with  a  numerical  example. 


4.  The  Survival  Probability  Function. 

In  the  present  section  we  establish  the  integral  equation  (1.1).  Let  Px  and  Pj,  be  a 
visible  point  on  C  and  a  point  to  its  right,  L*  <  x  <  y  <  U* .  The  Hunter  starts  shooting 
trials  when  the  target  is  at  Px .  The  attack  terminates  when  the  target  reaches  Py ,  if 
it  has  not  been  destroyed  before.  Let  S(x,  y)  designate  the  survival  probability  function. 
We  recognize  three  exclusive  and  exhaustive  events. 

(i)  The  visible  segment  to  the  r.h.s.  of  Pr  terminates  at  a  point  to  the  right  of  Py ; 

(ii)  The  visible  segment  on  the  r.h.s.  of  Px  terminates  at  a  point  P<,  t  <  y,  and  the 
length  of  the  shadow  starting  at  Pt  is  greater  than  y  —  t. 

(iii)  The  visible  segment  on  the  r.h.s.  of  P*  terminates  at  a  point  Pt,  t  <  y,  and  the 
length  of  the  shadow  to  the  r.h.s.  of  P*  is  smaller  than  y  —  t. 

As  mentioned  in  Section  1,  the  conditional  survival  probability  of  a  target  moving  on  a 
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visible  segment  of  length  L  is  exp {—q,L},  for  some  0  <  q  <  oo.  Accordingly, 
5(x,y)  =  C-^-*>(l-K(y-x|x)) 


(4.i; 


+  J*  -  D(y  |  t))dV(t  -  x  \  x) 


+ 


d 


J 


f  S(z,y)D'(z\t)dz\dV(t-x\x), 


where  D'(z  |  t)  =  —D(z  \  t)  is  the  p.d.f.  of  D(z  \  t).  Notice  that  D'(z  |  t)  =  0  for  all 
t  <  z  <  um(t )  +  t.  Let  zm(t)  =  um(t)  +  t.  zm(t )  is  the  first  term  on  the  r.h.s.  of  (3. IS) 
with  x  =  t.  Let  tm(z)  be  the  inverse  of  zm(t )  then,  by  changing  the  order  of  integration 
we  obtain 


(4.2 


j: 


,-q(t-x) 


X  |  x) 


=j: 


|  r  S(z,y)D'(z  |  t)dz  1  dV(t  — 

Utm(z)-X  ) 

e  qtD'(z  I  t  +  x)dV(t  I  x)  >  dz. 

Accordingly,  define 

(4.3)  A(x,  y)  =  e-^y~x\l  -V(y-x\  x))  +  *  e~qt(l  -  D{y  \  t  +  x))dV{t  \  x), 

Jo 

and 

r^m{z)  —  x) 

(4.4)  B(x,z)=  /  e~qtD'(z  (  t  +  x)dV{t  \  x). 

Jo 

Thus,  the  integral  equation  (4.1)  can  be  written  as  in  (1.1). 


5.  Algorithms  For  Discrete  Approximations  And  Numerical  Examples. 

In  the  present  section  we  consider  discrete  approximations  to  the  functions  Hn(t\x), 
Gn(ti,t2',x),  n  >2  and  S(x,y). 

For  a  given  integer,  N ,  partition  the  interval  (x,y)  to  N  subintervals.  Accordingly,  let 
6  —  (y  -  x)/N ,  t0  =  x  and  tj  =  t0  +  jf>,  j  =  0,  !,•••  ,N. 
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For  i  =  0,  ■  •  •  ,  -V ,  let 


(5.1) 


Hi(i)  =  H(ti;  t0)  —  exp{~n(to,  *,)}• 


We  compute  next  the  function  i),  =  0,---  ,N ;  which  is  a  discrete  approximation 

to  (3.13).  For  i  =  0,  -  ■■  ,N  and  j  =  i,  ■  ■  •  ,  N , 

i 

(5.2)  =  ^exp ■  [Hx(k)  -  Hx{k  -  1)], 

fc=0 

where  Hx(  —  1)  =  0.  Notice  that  -  Hx(0)  for  all  j  =  0,1,---  ,N;  and  for  i  >  1, 

j  >  i 

i 

(5.3)  G2(i,j)  =  #a(0)  +  expM*o,t,)}  •  ^ exp{-M<*>  tj)){Hi(k)  -  Hi(k  -  1)]. 

Jt=i 

Moreover, 


(5.4)  G2(i,j)  =  G2(j,j)  for  all  i  >  j. 

We  compute  afterwards  recursively,  for  every  n  >  3 ,  i  =  -  ,1V,  j  =  i,  -  ■  •  -IV 

Gn(i,j)  =  y^y^exp{-(/i(t/,t>)  - -  Gn-i(k  -  1,1) 

(5-5)  fc=0  /=fc 

—  Gn-l(k ,  1—1)4-  Gn-l(k  —  1, 1  —  l)]j 

and  for  i  >  j 

(5.6)  Gn{hj)  ~  Gn(j\  j)- 

For  1  =  0,  Gn(0,j)  =  Gn_i (0,  j )  =  tfi(0),  j  =  0,---  ,1V.  After  computing  these  functions 
we  determine  =  G»Cm),  J  =  0,1,-*-  ,1V.  Fn(j)  is  the  discrete  approximation  to 

the  c.d.f.  of  Tn(x),  namely  Hn(t;x );  i.e.,  Hn(tj-,x)  «  Hn{j). 

In  Table  5.1  we  present  numerical  results  obtained  by  applying  this  algorithm  to  the 

following  special  case. 
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We  consider  a  standard-uniform  Poisson  field,  with  uniform  distribution  for  the  disk 
radius  on  the  interval  (a,  b) .  In  the  appendix  we  present  the  functions  K±(x,t),  t  >  0,  for 
this  case.  We  compute  the  numerical  example  for  Table  5.1  with  the  following  geometrical 
parameters:  r  —  100[m],  u  =  40[m],  w  =  60[m],  a  =  l[m],  b  =  2.5[m],  x  =  10[m], 
I*  =  —  100[m],  U*  =  100[m],  We  present  in  the  tables  the  values  of  Hn(j ),  n  =  1,2.3, 
j  —  0,  •  ■  ■  ,20,  when  6  =  l[m] . 


Table  5.1.  Values  of  Hn(j)  for  two  values  of  A. 


A  =  0.02  [1/m2] 

A  =  0.2  [1/m2] 

j 

HxU) 

H2(j) 

Hs(j) 

HiU) 

H2(j) 

Hs(j) 

0 

0.4914 

0.4914 

0.4914 

0.0000 

0.0000 

0.0000 

1 

0.5434 

0.5434 

0.5434 

0.0000 

0.0000 

0.0000 

o 

0.6009 

0.600S 

0.6008 

0.0000 

0.0000 

0.0000 

3 

0.6645 

0.6642 

0.6442 

0.0003 

0.0003 

0.0003 

4 

0.7342 

0.7336 

0.7336 

0.0021 

0.0020 

0.0020 

5 

0.S039 

0.S030 

0.8030 

0.0130 

0.0119 

0.0119 

6 

0.S664 

0.8652 

0.8652 

0.0576 

0.0522 

0.0522 

7 

0.9182 

0.9170 

0.9170 

0.1828 

0.1663 

0.1663 

8 

0.9569 

0.9559 

0.9559 

0.4164 

0.3843 

0.3842 

9 

0.9811 

0.9804 

0.9804 

0.6842 

0.6484 

0.6482 

10 

0.9934 

0.9930 

0.9930 

0.8763 

0.8507 

0.8505 

11 

0.9985 

0.9984 

0.9984 

0.9707 

0.9590 

0.9588 

12 

0.9999 

0.9999 

0.9999 

0.9981 

0.9958 

0.9957 

13 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

14 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

15 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

As  seen  in  Table  5.1,  the  convergence  of  Hn{j )  to  the  c.d.f.  of  U(x)  is  quite  rapid.  We 
have  therefore  approximated  the  c.d.f.  D(u  j  x)  by  the  sequence  D(j  j  i)  —  D(tj  |  t,-), 
i  =  0, 1,  •  •  •  ,  N ,  j  =  i,  i  +  1,  •  •  •  .  The  function  A(x,  y)  was  computed  for  the  arguments 
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tt ,  tj .  by  the  approximation 


i(JV,  JV)  =  1 

i(.V  -  1,  JV)  =  -  V(\  |  iV  -  1))  +  (1  -  D(N  |  iV  -  1)) 

(5.7)  M*  -  j,  N)  -  -  V(i  !  JV  -  ;)) 

i 

+  Z  j  JV  -  j)  -  V(f  -  1  |  JV  -  j))- 

i=i 

•  (1  —  D(N  \  N  —  j  +  /)),  j  =  2,  •  •  •  ,  JV 

where 

(5.8)  D(iV  |  JV  -  i)  =  l(D(N  |  N~i)  +  D(N  \N-i  +  1)), 


for  all  i  =  1,2.  •  •  •  ,  'V.  Recall  that  D(N  |  N)  =  0  =  D(N  |  JV  +  1). 
Similarly,  we  define 


(5.9) 


B(N,  JV)  =  0 

B(N  -  1,JV)  =  Ie-9i/2F(l  |  JV  -  1)D(JV  |  N  -  I), 


and  for  j  =  2,  •  •  •  ,  JV ,  /  =  1,  *  •  •  ,  j  we  compute 


(5.10) 


/ 

B(N  -  j,N  -j  +  l)  =  Yl e_?6(i"  |  iV  -  j)  -  V(i  -  1  |  JV  -  ;)]• 

i=  1 


•  [D(N  -  j  +  / 1  JV  -  j  +  i)  -  D(N  -  j  +  l-  l\  N-  j  +  *)], 


where  D(N  —  j  +  l  —  1  |  JV  —  j  +  /)  =  0. 

Using  these  sequences,  we  compute 

5(JV,  JV)  =  1 

5(JV  -  1,  JV)  =  A(N  -  1,  JV) 

and.  for  j  —  2,  ■  •  •  ,  JV 

N-l 

(5.11)  5(JV  -  j,  JV)  =  A(JV  -  j,  JV)  +  B(JV  -  j,  JV)  +  Y,  B(n  -  j,i)S(i,N). 

i=N —j+l 
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The  function  S(x,y)  is  approximated  by  5(0, N).  In  Table  5.2  we  present  the  values  of 
5( -V  —  j,  N) ,  for  the  geometrical  parameters  of  Table  5.1,  with  a  =  2[m] ,  b  =  3.5[m]  and 
several  values  of  A.  Also  here  6  =  l[m].  The  value  of  q  is  —ln( 0.8).  This  corresponds  to 
the  situation  in  which  one  shooting  trial  takes  as  long  as  the  target  travels  l[m],  and  the 
probability  of  destroying  the  target  in  one  trial  is  0.8. 


Table  5.2.  Survival  Probabilities  S(N  —  j,N ),  for  A  =  0.02(0.01)0.05. 


N  =  20;  a  =  2,  b  =  3.5,  r  =  100,  u  =  40,  w  =  60. 


AA 

0.02 

0.03 

0.04 

0.05 

0 

1.00000 

1.00000 

1.00000 

1.00000 

i 

0.81711 

0.82447 

0.83112 

0.S3715 

2 

0.69730 

0.72040 

0.74052 

0.75808 

3 

0.61877 

0.65867 

0.69187 

0.71965 

4 

0.56729 

0.62203 

0.66572 

0.70094 

5 

0.53352 

0.60025 

0.65163 

0.69179 

6 

0.51135 

0.58728 

0.64401 

0.68729 

7 

0.49677 

0.57954 

0.63987 

0.68505 

8 

0.48718 

0.57489 

0.63759 

0.68390 

9 

0.47460 

0.56551 

0.63027 

0.67817 

10 

0.45820 

0.55111 

0.61788 

0.66784 

11 

0.43776 

0.53154 

0.60001 

0.65204 

12 

0.41376 

0.50726 

0.57671 

0.63029 

13 

0.38732 

0.47959 

0.54917 

0.60348 

14 

0.36009 

0.45064 

0.51979 

0.57418 

15 

0.33382 

0.42283 

0.49155 

0.54589 

16 

0.30991 

0.39803 

0.46678 

0.52137 

17 

0.28902 

0.37709 

0.44652 

0.50189 

18 

0.27122 

0.35999 

0.43073 

0.48741 

19 

0.25614 

0.34607 

0.41851 

0.47687 

20 

0.24309 

0.33427 

0.40849 

0.46864 

13 
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Appendix 

The  Functions  K±(x,t)  in  the  Standard- Uniform  Case,  With  Uniform  Distribution  of 
radii  on  (a,  b) . 

Let  K+(x,  t ,  z )  denote  the  area  of  the  set  bounded  by  the  line  Ct ,  the  ray  Rr+t ,  t  >  0, 
and  the  lines  U  and  W;  £+  is  the  line  parallel  to  Rx ,  on  its  r.h.s.,  of  distance  2  from 
it.  This  is  the  set  of  all  disk  centers  between  Rr  and  Rr+t,  of  radius  Z  =  z,  which 
do  not  intersect  Rr.  In  order  to  simplify  notation,  we  assume  that  w  —  r.  In  actual 
computations  we  substitute  xw/r  and  tw/r  for  x  and  t  in  the  formulae  given  below.  Let 
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d  =  (x-  -f  tc2)1/2  .  Simple  geometrical  considerations  yield: 


(-4.1) 


\  _  f  zd  \  \w2  —  u2  zd 

R  +  {x,t,z)  =  1 1  —  <u\  — - t-  2— — 

{  t  J  [  2w  w  +  u 


,  zd 

+  I{u<  —  <w 


2tw 


( tw  —  zd)2 


where  I  {A}  is  the  indicator  set  function,  which  assumes  the  value  1  if  A  is  true,  and  the 
value  0  otherwise. 

Notice  that  A'+(x ,  t,  z)  depends  on  x  only  via  x2  .  Symmetry  implies  that  A'_(  —  x,  t ,  z)  = 
K+(x,t,z)  =  K+(—x,t-z)  for  all  —  cc  <  x  <  oo.  Hence,  K+(x,t)  =  K-(x,t)  and  we 
delete  the  ±  subscript  of  K.  Finally,  K(x,  t)  =  E{I\(x,t,  Z)}  with  respect  to  the  uni¬ 
form  distribution  of  Z  over  ( a.b ).  Let  xi  =  tu/d  and  xi  =  tw/d.  The  function  I\(x.t) 
assumes  the  following  forms: 

(i)  If  b  <  xi , 


w2  —  u2  f  d 

(.4.2)  /COM)  -  ( t  -  -M_(a  +  6 ) 

2u;  \  u  +  w 


(ii)  If  a  <  Xi  <  b  <  X2 


K(x,  t)  = 


2  2 
W  —  U 


(.4.3) 


+ 


2w 

1 


•T  j  CL 


b  —  a  u  +  w  b  —  a 


(*?  -  «2) 


d2 


2  2^  ^  / 12  2\  , 
t*w  - - tu>- - (6  —  xt)4- 

2tw  V  b  —  a  b  —  a  3(6  —  a) 


(63-*2) 


(iii)  If  a  <  Xi  <  X2  <  6 


(-4.4) 


+ 


2u; 

I 


6  —  a  u  +  w  b  —  a 


(xi  —  a2) 


2tw 

(iv)  If  xx  <  a  <  6  <  x2 , 


1  -tw-r^—{x l  -  x\)  + -—^—-{xl  -  x\) 

I)  —  a  6  —  a  o(o  —  a) 


2  2 12  "^1 

t  w 
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(-4.5) 


tw  d  .  d?(a?  +  ab  +  b2) 

A(l-‘>=T-2(a  +  6)+  Cta - - 


(v)  If  xi  <  a  <  X2  <  b 


(.4.6)  A-(x. «)  =  -  ^(xi  -  ^  +  ^^(4  -  A 

(vi)  If  X2  <  a 


(.4.7) 


I\(x,  t )  =  0. 


16 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Dmtm  En fr.d) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER  2.  GOVT  ACCESSION  NO. 

A-**  25W  ./-M 

3.  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  (end  Subtitle) 

THE  SURVIVAL  PROBABILITY  FUNCTION  OF  A  TARGET 
MOVING  ALONG  A  STRAIGHT  LINE  IN  A  RANDOM  FIELD 

OF  OBSCURING  ELEMENTS 

5.  TYPE  OF  REPORT  6  PERIOD  COVERED 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  author^.; 

S.  Zacks  and  M.  Yadin 

8.  CONTRACT  OR  GRANT  NUMBERC*; 

DAAL03-89-K-0129 

9-  PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 

State  University  of  New  York 

Binghamton,  NY  13901 

10.  PROGRAM  ELEMENT.  PROJECT.  T  ASK 
AREA  6  WORK  UNIT  NUMBERS 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

U.  S.  Army  Research  Office 

Post  Office  Box  12211 

Research  Triangle  Park.  Nr  ?77fN 

12.  REPORT  OATE 

May  15,  1990 

13.  NUMBER  OF  PAGES 

16 

14.  MONITORING  AGENCY  NAME  ft  ADDRESS (II  dlllerent  from  Controlling  Office) 

IS.  SECURITY  CLASS,  (ot  thle  report) 

Unclassified 

15a.  DECLASSI  FI  CATION/ DOWN  GRADING 
SCHEDULE 

16.  DISTRIBUTION  STATEMENT  (ol  thle  Report) 

Approvcb  lor  public  release;  «ributio» 

unlimited. 

17.  DISTRIBUTION  STATEMENT  (of  the  ebetrect  entered  In  Block  20,  II  dlllerent  from  Report) 

16.  SUPPLEMENTARY  notes 

•rum  ..'rW.  OPINIONS,  AND/OR  FINDINGS  CO'-”  rn  ^ 

-  *032  OF  THE  AUTHORS)  AND  SHOULD  *■'  AS 

/■F  CIAL  DEPARTMENT  OF  THE  ARMY  POOMiCM,  F— 'C  .  CH  UE 

CISION.  UNLESS  SO  DESIGNATED  BY  OTHER  DOCUMENTATION, 

19.  KEY  WORDS  (Continue  on  rereree  elde  If  neceeeery  end  Identity  by  block  number) 

Line  of  sight;  visibility  probability;  distribution  of  shadow  length; 
survival  probability;  Poisson  random  field. 

ZO.  AHTFACr  (CmrOmum  me  rereree  etde  M  neceeeery  end  Identity  by  block  number) 

A  target  is  moving  along  a  straight  line  path.  Random  portions  of  the  path  might 
je  invisible  to  the  hunter  (in  shadow).  Shooting  trials  axe  conducted  only  along  the 
risible  segments  of  the  path.  An  algorithm  for  the  numerical  determination  of  the  survival 
probability  of  the  target  is  developed.  This  algorithm  is  based  on  the  distribution  of 
ihadow  length  which  is  also  developed. 

i 

DO  ,  jSTn  U73  COITION  OF  t  NOV  «*  IS  OBSOLETE 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Dele  Bnleted) 


t 


